First determination of the strange and light quark masses from full lattice QCD 
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We compute the strange quark mass m s and the average of the u and d quark masses m using 
full lattice QCD with three dynamical quarks combined with the experimental values for the n and 
K masses. The simulations have degenerate u and d quarks with masses m u — irid = rh as low as 
m s /8, and two different values of the lattice spacing. The bare lattice quark masses obtained are 

converted to the MS scheme using perturbation theory at 0(as)- Our results are: m^ IS (2GeV) = 
76(0) (3) (7) (0) MeV, rh H5 (2GeV) = 2.8(0)(1)(3)(0) MeV, and m s /rh = 27.4(1)(4)(0)(1), where the 
errors are from statistics, simulation, perturbation theory, and electromagnetic effects, respectively. 



I. INTRODUCTION 



The masses of the strange and light quarks are fun- 
damental parameters of the Standard Model that are a 
priori unknown and must be determined from experi- 
ment. This is complicated, however, by confinement in 
QCD, so that quarks cannot be observed as isolated par- 
ticles. We can only determine their masses by solving 
QCD for observable quantities, such as hadron masses, 
as a function of the quark mass. This can be accom- 
plished with the numerical techniques of lattice QCD. 
Precise knowledge of quark masses constrains Beyond 
the Standard Model scenarios as well as providing in- 
put for phcnomcnological calculations of Standard Model 
physics. The strange quark mass, in particular, is needed 
for various phenomenological studies, including the im- 
portant CP-violating quantity e'/e where its uncer- 
tainty severely limits the theoretical precision. 

Previously, shortcomings in the formulation of QCD 
on the lattice and limitations in computing power have 
meant that lattice calculations were forced to work with 
an unrealistic QCD vacuum that either ignored dynam- 
ical (sea) quarks or included only u and d quarks with 
masses much heavier than in Nature. This condemned 
determinations of the quark masses to rather large sys- 
tematic errors (10-20%) arising from the inconsistency 
of comparing such a theory with experiment. The de- 
termination presented here uses simulations with the im- 
proved staggered quark formalism that have a much more 



realistic QCD vacuum with two light dynamical quarks 
and one strange dynamical quark. We describe how the 
bare quark masses in the lattice QCD Lagrangian can 
be fixed using chiral perturbation theory to extrapolate 
lattice results to the physical point, and how the lat- 
tice quark masses obtained can be transformed to a con- 
tinuum scheme (MS) using lattice perturbation theory. 
Working in the region of dynamical u/d quark masses 
below m s /2 and down to m s /8 gives us control of chi- 
ral extrapolations and avoids the large systematic errors 
from dynamical quark mass and unquenching effects that 
previous calculations have had. 

Staggered quarks are fast to simulate. They keep a 
remnant of chiral symmetry on the lattice, and therefore 
give a Goldstone pion mass which vanishes with the bare 
quark mass. This allows the relatively simple determi- 
nation of the quark mass described here, which is not 
available, for example, in the Wilson quark formalism. 

The staggered quark formalism does have several un- 
wanted features, however. With the naive staggered ac- 
tion, large discretization errors appear, although they are 
formally only 0(a 2 ) or higher (a is the lattice spacing). 
The renormalization of operators to match a continuum 
scheme can also be large and badly behaved in pertur- 
bation theory. This is true, for example, for the mass 
renormalization that is needed here. It turns out that 
both problems have the same source, a particular form of 
discretization error in the action, called "taste violation," 
and both are ameliorated by use of the improved stag- 
gered formalism 0] ■ The perturbation theory then shows 
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small rcnormalizations jiM M an d discretization errors 
are much reduced [(J 0, 0- Empirically, taste violation 
remains the most important discretization error in the 
improved theory, despite being subleading to "generic" 
discretization errors. The Goldstone meson masses we 
will discuss here are affected by this at one-loop in the 
chiral expansion. Staggered chiral perturbation theory 
(SXPT) [1 El [HI 03 allows us to control these effects 
and reduce discretization errors significantly. 

A more fundamental concern about staggered fcrmions 
is based on the need to take the fourth root of the 
quark determinant to convert the four-fold duplication 
of "tastes" into one quark flavor. It is possible that there 
are nonlocalities in the continuum limit that would spoil 
the description of QCD at some level. Checks of the for- 
malism against experimental results 0, IT^l IbH ITEL ITfil | , 
make this unlikely, we believe, but further work along 
these lines is crucial and continuing. 



II. LATTICE DATA 

The simulation data of the MILC collaboration 0,0] 
are analyzed' staggered quarks with leading errors at 
0(a s a 2 ,a 4 ) □ and one-loop Symanzik improved gluons 
with tadpole-improvement [nj- Two sets of con- 
figurations are used: a "coarse" set at lattice spacing 
a « 1/8 fm and sea quark masses of am' u = am' d = 
am' = 0.005, 0.007, 0.01, 0.02, 0.03 with am' s = 0.05, 
and a "fine" set at a w 1/11 fm with sea quark masses 
of am' = 0.0062, 0.0124 and am' s = 0.031. Here we use 
primes on the sea quark masses to emphasize that these 
are the nominal quark masses used in the simulation, not 
the physical masses m s or m = (m u + m^jl. The simu- 
lations are "partially quenched," with a range of valence 
masses from m' s down to m' s /10 (coarse) and m' s /5 (fine), 
not necessarily equal to the sea quark masses, simulated 
on each lattice. It should be noted that the quark masses 
in lattice units quoted here contain a factor of mop, the 
tadpole-improvement factor determined from the fourth 
root of the average plaquette, compared with a more con- 
ventional definition of quark mass 2] . This is taken care 
of nonperturbatively before our renormalization below. 

The lattice spacing a is determined ultimately from 
the T'-T mass difference a useful quantity because 
it is approximately independent of quark masses, includ- 
ing the 6-mass. An analysis of a wide range of other 
"gold-plated" hadron masses and decay constants on 
these configurations shows agreement with experiment at 
the 2-3% level ^3|. Gold-plated hadrons are stable (in 
QCD), with masses at least 100 MeV below decay thresh- 
olds, so their masses are well-defined both experimentally 
and theoretically, important for fixing the parameters of 
QCD. The only gold-plated light mesons available to fix 
m and m s are the tt and K. There is none with only s 
valence quarks because the <f> is unstable and the pseu- 
doscalar is strongly mixed. Baryons can provide an al- 
ternative, the nucleon for m and the Q for m s , but their 



statistical errors are large, and they are not very sensitive 
to the quark masses. 

Our analysis uses SXPT ^] to fit the dependence of 
the results on the quark masses. This dependence can 
then be extrapolated/interpolated to the point where the 
(Goldstone) tt and K have their physical masses, thereby 
determining the bare lattice m and m s . At the level of 
precision at which we are working, and because we take 
m u = md-, we must be careful about electromagnetic 
(EM) and isospin- violating effects. At lowest nontrivial 
order in e 2 and the quark masses, Dashen's theorem [2lJ 
states that m? + and m 2 K+ receive equal EM contribu- 
tions; while the tt° and K° masses are unaffected. How- 
ever, at next order, there can be large and different con- 
tributions to and m 2 K+ of order e 2 m 2 K [22l I2H. l24l] . 
Let Ae parameterize violations of Dashen's theo- 
rem: (m 2 v + - m 2 )em = (l + A_ E )(m 2 r+ -rn^EM- Then 
Refs. milEl suggest A E w 1. 

Including EM and isospin effects, the physical values 
of ?7i and m s can then be determined by extrapolating 
the lattice squared meson masses to m| = m 2 and 
m \ = ( m K° + m K+ - C 1 + A E)(ml + - m 2 o))/2, us- 
ing experimental values on the right hand side of these 
expressions. We are neglect ing O (( m u — md) 2 ) correc- 
tions, which should be tiny [2(j. EM contributions to 
the neutral particle masses are also neglected, and we 
take account of this in our error. For the tt° the violation 
of Dashen's theorem is 0(e 2 m 2 /(8ir 2 f%)) and negligible. 
For m 2 K0 the violation is in principle the same order as 
for m 2 K+ [23, but in model calculations [24[ it appears 
to be very small. To be conservative, we consider EM 
contributions to m 2 K0 of order of half the violations of 
Dashen's theorem, with unknown sign. Effectively, this 
replaces Ae w 1 in the formula for m 2 ^ above with the 
range 0-2, which we take as the EM systematic error. 

III. CHIRAL FITS AND SYSTEMATIC ERRORS 

Here we briefly describe the fits to SXPT theory forms 
and the estimate of the associated errors EEl. Be- 
cause the squared meson masses (M^eson) are nearly lin- 
ear in the valence quark masses, the final values of the 
quark masses are quite insensitive to details of the chiral 
fits. Chiral logs and NLO (and higher) analytic terms 
only affect the results at the s=s 5% level. 

SXPT is a joint expansion in x q and x a i , which arc 
dimensionlcss measures of the size of quark mass and 
lattice spacing effects, respectively: 

_ 2[im q a 2 A 

Xq ~ s^/f ; Xa2 = 8^7| ' ( ' 

m q is the quark mass, 1[im q is the tree-level mass of a 
qq meson, and / w f=a 131 MeV. a 2 A is an average meson 
splitting between different tastes. On the coarse lattices 
x a 2 « 0.09; on the fine, x a i « 0.03. 

For physical kaons, the relevant expansion parameter is 
x u d tS = (x u d+x s )/2 « 0.18. Since our lattice data is very 
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precise (0.1 to 0.7% on M^ cson ), it is clear that we cannot 
expect NLO or even NNLO XPT to work well up to the 
kaon mass. If however the valence quark masses are lim- 
ited by m x + m y < 0.75m^, we obtain good fits including 
NNLO analytic terms. Such fits are consistent with XPT 
expectations: the coefficients of NLO and NNLO terms 
are 0(1) when these terms are expressed as functions of 
x q and x a 2 . When fitting up to the strange mass we in- 
clude NNNLO as well as NNLO terms, but satisfy the 
chiral constraints by fixing the NLO terms from lower 
mass fits. Since the s quark mass can be reached in sim- 
ulations, the form of the NNLO and NNNLO terms is 
not important; such terms simply allow for a reasonable 
interpolation to the physical m s . 

Both decay constant and M^, eson data and both coarse 
and fine ensembles are fit simultaneously. Although 
NLO taste-violations are explicitly included, we allow 
for "generic" discretization errors by using a Bayesian 
fit [27| that permits physical parameters to change by 
order asa 2 Aq CD ~ 2% in going from the coarse to the 
fine configurations. 

The T system provides an absolute lattice scale, but it 
is convenient to use the relative scale determined from 
n, a parameter derived from the heavy quark poten- 
tial [23, , to compare accurately the scale for different 
sea quark masses within the coarse or fine set. T split- 
tings give ri = 0.317(7) (3) fm [l4|. Using the volume de- 
pendence calculated in NLO SXPT Jl3ll3, ( and tested 
against results on different volumes [lj|) the small finite- 
volume effects (< 0.75% in M^ eson ) can be removed from 
our data with negligible residual error. 

Figure ^ compares our fit with our partially quenched 
data for M^ eson . The data appear quite linear to the 
eye. Indeed, linear fits change our result for the quark 
masses by only 2 to 7%, depending on the fit range chosen 
and whether or not the correlated decay constants are fit 
simultaneously. However, since the statistical errors in 
our data are so small, the nonlincarities from chiral logs 
and higher order analytic terms are crucial for obtaining 
good fits: linear fits have x 2 /dof~ 20. Nonlinear fits have 
a confidence level of 0.28, are crucial to obtaining Gasser- 
Leutwyler parameters and affect the decay constants by 
-4-12%. 

We extrapolate/interpolate in mass on the coarse and 
fine lattices separately to find the lattice values of the 
light and strange masses that give m| and m^.. We 
get am s = 0.0390(1)(20), am = 0.00141(1)(8) on the 
coarse lattices and 0.0272(1)(12) and 0.000989(3)(40) on 
the fine, where errors are statistical and systematic. The 
systematic errors are dominated by the chiral extrapo- 
lation/interpolation, estimated by varying the fits, and 
the scale uncertainty (EM effects account for the slight 
difference with 0|). Alternatively one can extrapolate 
the chiral fit parameters to the continuum, setting taste- 
violating parameters zero, and then perform the chiral 
extrapolation/interpolation to the physical masses. This 
is shown as the dashed green lines in Fig. ^ The meth- 
ods give final MS masses that differ by less than 2%. We 
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FIG. 1: Partially quenched data for squared meson masses 
made out of valence quarks x and y as a function of m x /m' s . 
We show results from two lattices: a coarse lattice with sea 
quark masses am' — 0.01, am' s — 0.05, and a fine lattice with 
am' = 0.0062, am' s = 0.031. Three sets of "kaon" points with 
m y = m' s , 0.8m' s , 0.6m' s , are plotted for each lattice. "Pion" 
points have m x = m y . The solid lines come from a fit to all 
the data (not just that plotted). The statistical errors in the 
points, as well as the variation in the data with sea quark 
masses are not visible on this scale. The green dashed lines 
give the continuum fit described in the text, and the magenta 
vertical dotted line gives the physical m/m s obtained. 

choose the first method for the central values and include 
the variation with method in the systematic error. 

The same SXPT fits that produce the quark masses 
above give Gasser-Leutwyler parameters in reasonable 
agreement with phenomenological valu es |l2l and j v and 
Jk in agreement with experiment |l2t Il3l|. Final results 
and all details of the fits will be described in Ref. 0| . 

It is important to provide further checks of m s and m 
using other gold-plated masses and mass differences. We 
focus on m s because it has smaller statistical error and 
less dependence on chiral extrapolations. From the heavy 
hadron sector 2ms aqj , s — mx is sensitive to m s but not to 
other masses. Here 2mB av . s is the B s , B* spin-averaged 
mass, used to reduce dependence on the coefficients of 
relativistic corrections in the 6-quark action. Note, how- 
ever, that the B* is close to decay threshold and may not 
be gold-plated. Figure|2]shows coarse-lattice data for this 
splitting. The results are 2% high, but this is also our 
estimate of discretisation errors in the calculation (we do 
not expect sensitivity to taste-violation |30jh This quan- 
tity then provides a check of our m s determination at the 
20% level because the experimental splitting varies only 
by w 15% in changing from m to m s . Figure[21also shows 
results for the £1 baryon mass, on both coarse and fine 
ensembles. Although statistical errors are large there is 
a trend downwards on the finer lattices and signs that 
a continuum extrapolated result will agree with experi- 
ment. An expected 2% error on the final value for tuq 
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FIG. 2: Lattice results for two masses which show sensitivity 
to m a , plotted against rh' /m' s . The valence s masses are taken 
at the m s values determined here. The bursts give the corre- 
sponding experimental result. The squares are 2m_B s „„ — mr 
for two of the coarse ensembles. The upper results are for the 
mass of the fl (sss) baryon, on both coarse (diamonds) and 
fine (crosses) ensembles. 



would lead to a 6% determination of : 



where b(am) w 0.5432 — OAQ(am) 2 , correct to 0.1% up 
to (am) = 0.1. 7o = — is the universal one-loop anoma- 
lous mass dimension. Naive staggered quarks have a 
poorly convergent Z m with 6(0) ~ 3.6 as a result of 
taste- violations. It is clear that the improved staggered 
quark result is much better. Tadpole-improvement is 
also important, because of the long paths of gluon fields 
required to suppress taste-violations. Without tadpole- 
improvement 6(0) = 2.27. 

We match our lattice to the MS scheme at the target 
scale of 2 GeV, though the results and errors are not sen- 
sitive to this choice. Because the mass renormalization 
has an anomalous dimension, the optimal q* value for ay 
at this scale is dependent on a. q* is set by a second order 
BLM method On the fine lattices, q* is 1.80/a HJ 
and av(q*) = 0.247(4) in Z m . On the coarse lattices, 
q* = 2.335/a, giving ot\r{q*) = 0.252(5). A conservative 
estimate of the perturbative error in Z m , informed by the 
chiral fits, is 1.5 x a 2 , «9%. 

This gives values of 74.3 MeV on the fine lattices 
and 72.3 MeV on the coarse lattices. Our central values 
are obtained by extrapolating linearly in ago 2 , the size of 
the leading discretization errors. Alternative's, such as a 
linear extrapolation in a s a 2 , the size of taste-violations, 
or a continuum-extrapolated chiral fit, give results that 
vary by less than 1 MeV, which we take as the extrapo- 
lation error and fold into the total systematic error. Our 
final quark masses are: 



MS 



IV. CONNECTING m lattice WITH m MS 

The continuum quark mass in the conventional modi- 
fied Minimal Subtraction scheme is determined from: 

m Ms (//)= (^)o (i +av{q * )Z (2) {atjLjiam)o) + o(a 2 )) , 

(2) 

where (am)o is the a posteriori tuned bare mass in lattice 
units obtained above, converted from the MILC conven- 
tion by dividing by itop. Z m is the mass renormalization 
that connects the bare lattice mass and the MS mass. 
The strong coupling constant in the V scheme is set using 
third order perturbative expressions for the logarithms of 
small Wilson loops 0, [3^ compared with lattice results 
on these configurations. The value obtained is run to an 
optimal scale q*, chosen as described below. 

Z m is calculated by connecting the bare quark-mass 
to the pole-mass in lattice perturbation theory and 
using the pole mass to MS mass relation |3^| at one loop. 
The lattice calculation was done both by hand and us- 
ing automated methods 0, H^, which become increas- 
ingly important for improved actions. The evaluation has 
been checked to lower precision via a completely differ- 
ent method [3(| . Integrals were evaluated here using the 
numerical integration package, VEGAS [13 ■ We find 

Z$(an, am ) = (b(am ) - - ln(a/x) ] , (3) 
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(2GeV) = 76(0) (3) (7) (0) MeV (4) 

(2GeV) = 2.8(0) (1) (3) (0) MeV (5) 
m s /m = 27.4(1) (4) (0)(1) , (6) 



where the errors come from statistics, simulation system- 
atics, perturbation theory, and electromagnetic effects, 
respectively. The systematic error includes the scale er- 
ror in quadrature with the chiral and continuum extrap- 
olation errors. The ratio m s /rh in Eq. @ is almost in- 
dependent of the perturbation theory. It is also strongly 
constrained by the fact that 2m 2 K — m 2 is almost in- 
dependent of light quark mass over a large range. For 
our coarse lattices it increases by 2% as m' changes from 
m' s /5 to m' s ; for the fine lattices by 4%. 



V. COMPARISON WITH PREVIOUS 
DETERMINATIONS 

There is a long history of sum rule determinations of 
the strange quark mass, with the gene ral trend of de- 
creasing values. The current status [33, l40l l4lj is broad 
agreement between results from scalar and pscudoscalar 
spectral functions and from SU(3) breaking in t hadronic 
decays, with m s around 100(20) MeV. The latter method, 
however, is sensitive to the value of |V^ S |. Lattice results 
in the quenched approximation give values around 100 
MeV but more recent results with two flavors of rather 
heavy dynamical quarks give a smaller value around 90 
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MeV 01 ■ Both quenched and nj = 2 results suffer from 
the inherent systematic error of comparing an unphys- 
ical theory with experiment: results depend on what 
hadronic masses are used. Some determinations also do 
not use gold-plated quantities. JLQCD quote a pre- 
liminary rif = 3 result of 75.6(3.4) MeV, not yet including 
discretization and finite volume errors. They use clover 
quarks with rh! > m s /2, setting a with the p mass. 

Here we give results from nj — 3 simulations in the chi- 
ral regime. Using gold-plated quantities to fix the QCD 
parameters means that there is no remaining ambiguity 
in the match between QCD and experiment. The value 
we obtain for m s is lower than previous results, but we 
maintain that it is based on a firmer footing. It violates 
some quoted bounds from sum rules |44j . but these are 
open to question |4l| . Our result for m s /m is signifi- 
cantly larger than that determined from NLO XPT phe- 
nomenology , but is compatible with a NNLO analy- 
sis Gjs. We believe that existing staggered-quark results 
I12LI13LIT4IIT5I I make it unlikely that there are fundamen- 
tal problems with the formalism we are using. 
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